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ABSTRACT 



Linear state signal feedback is used to obtain exponential 
response from fourth order systems* Characteristic equation roots 
are selected to provide the desired exponential response with con- 
straint on initial conditions and system acceleration * A digital 
computer root locus program is developed to determine feedback 
coefficients in a manner which minimizes the possibility of oscilla- 
tory response in the presence of state sensor errors* The effect 
of noise on the state signal is investigated and a sample data fil- 
tering technique developed. A quasi-optimum time technique utilizing 
second order switching logic for initial control effort and linear 
feedback for terminal control is developed* 



ii 



ACKNOWLEDGMENT 



The author wishes to express his appreciation for the assistance 
and encouragement given him by Dr* H. A, Titus^ Jr oi> of the United 
States Naval Postgraduate School, 



iii 



I 



TABLE OF CONTENTS 



Section 




Title 




Pag 


1.0 




Introduction 




i 


2.0 




Linear Feedback for Exponential Settling 


2 




2.1 


Use of the Uncoupled Form 




4 




2.2 


Use of a Root Locus Study 




8 




2.3 


Consideration of Acceleration 


Constraint 


9 




2.4 


Statistical Considerations 




12 


3.0 




Use of Second Order Switching 
Quasi-Optimal Response 


Logic for 


18 


4.0 




Conclusions 




29 






Bibliography 




30 






APPENDICES 






I. 


Root Locus Study 




32 


II. 


Acceleration Constraint 




39 


III. 


Gaussian Noise Generator 




42 


IV. 


Statistical Study 




45 


V. 


Quasi-Optimal Control 




51 



iv 



LIST OF ILLUSTRATIONS 



Page 



Figure 



1. 


Block Diagram of Aircraft Longitudinal Motion 
with State Variable Feedback 


3 


2 . 


Block Diagram of Uncoupled Form of Aircraft 
Longitudinal Motion with State Variable 
Feedback 


7 


3. 


Block Diagram of General Fourth Order System 
with Initial Conditions on the State Variables 


10 


4. 


Acceleration of Fourth Order System with 
Dominant Characteristic Equation Root of 0,2 


14 


5. 


Time Response of Fourth Order System with 
Noisy Sensors 


15 


6. 


Time Response of Fourth Order System with 
Noisy Sensors 


17 


7 . 


Time Response for Saturating Linear Control 


23 


8. 


Time Response when Second Order Switching Logic 
is Applied to the Averaged Uncoupled States 


24 


9. 


Time Response when Second Order Switching Logic 
is Applied to the Uncoupled State Fair Having 
the Higher Velocity 


25 


10. 


Time Response when Second Order Switching Logic 
is Applied to the Uncoupled State Fair Having 
the Higher Energy 


26 


11 . 


Time Response Using Energy Priority and Revised 
Second Order Switching Logic 


27 


12. 


Time Response Using Energy Priority and Final 
Second Order Switching Logic 


28 



v 



1.0 Introduction 



The purpose of this study was to investigate means of controlling 
fourth order systems to obtain various types of responses. The initial 
portion of the study deals with the use of linear state signal feedback 
to obtain exponential settling. Use of the uncoupled form is investi- 
gated for use in the feedback solution analysis. The effect of noise 
on the state signals for the linear feedback solution was investigated 
and a sample data filtering technique developed. A root locus program 
was developed to determine feedback coefficients in a manner which will 
minimize the possibility of oscillatory response in the presence of 
state sensor errors. 

The remainder of the study deals with a quasi-optimum time 
solution which incorporates switching logic. Full control effort is 
used initially and switching takes place at predetermined levels of 
state variable combinations. After switching, control is allowed to 
decay for a linear termination of the solution. The objective here 
is to obtain a near- optimum time response with no possibility of con- 
trol chatter. 

Section 2.0 of this report contains the linear feedback portion 
of the investigation. Section 3«0 contains the quasi-optimum investi- 
gation. All synthesis for the study was conducted on the school 0 s 
Control Data Corporation 160k digital computer using fortran. All 
programming will be referenced in the text and shown in the appendices. 
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2.0 Linear Feedback for Exponential Settling 



The response of a system to a set of initial conditions may be con* 
trolled by feedback of the system state variables. That is, the forcing 
function of the system is made up of predetermined amounts of each of 
the system 1 s state variables. Thus, the characteristic equation of the 
controlled system may be adjusted to obtain the desired response. 

If an exponential response of the system's position state variable 
is desired, the characteristic equation will have only negative real 
roots, with a dominant root that causes the desired exponential path 
after the decay of the non-dominant roots. The proper feedback co- 
efficients of each of the state variables may be calculated to give the 
desired characteristic equation. 

Example 

An aircraft landing flare is representative of a class of auto- 
matic control problems in which a system has initial conditions of each 
state variable and it is desired that the position state variable 
settle to zero in an exponential manner. In order to design a suit- 
able control system, a mathematical description of the aircraft longi- 
tudinal motion is required. Assuming constant air speed and a shallow 
glide angle leads to the short period equation of longitudinal motion 



relating elevator position, d (radians), to aircraft pitch angle, 
p (radians ); 




These are written in terms of the following transfer function 




d(s) 




2 



where 



K = short period gain 



w = short period resonant frequency 
a = short period damping factor 
T = path time constant 

In order to complete the mathematical description of the aircraft, 
the transfer function relating altitude, h (feet), and pitch angle, 
p (radians), in terms of velocity V (feet per second) and path time 
constant T is 



h(s) = 



P(s) 



s(Ts + 1) 

Combining (1) and (2) results in a transfer function relating 
altitude and elevator position: 



( 2 ) 



h(s) = 



KV 



5^(s^ + 2as + w^) 



d(s) 



( 3 ) 



Letting the system forcing function u be equal to KVd^ a signal 
flow diagram of the system is shown in Figure 1. 




Figure 1. Block Diagram of Aircraft Longitudinal Motion with 
State Variable Feedback 
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The characteristic equation for the system of Figure 1 is 

s k + (b^ + 2a)s^ + (b^ + w 2 )s 2 + b>s + b^ ~ 0 {h) 

It is now possible to solve for the feedback coefficients 

bf, b^, by and b^ such that the characteristic equation will have 
the desired roots . Suppose, for a particular aircraft 
a = 0*5 and w =3 1.0, and the desired closed loop characteristic equa- 
tion roots are s =-0.l8, -1„0, -1.0, and -5.0s giving the charac- 
teristic equation 

s k + 7.18s 3 + 12 . 26 s 2 + 6 . 98 s + 0.9 » 0. ( 5 ) 

Equating (U) and ( 5 ) will then give the desired feedback co- 
efficients 

b L = 0 . 9 , b 2 = 6 . 98 , b 3 . 11 . 26 , b^ » 6.18. (6) 

This will give a time response solution for the flare, after allowing 

a short decay time (t^ ± 6 seconds) for the non-dominant roots, of 
approximately 

h(t) = h(t x ) e'°* l8t (7) 



2 . 1 Use of the Uncoupled Form 

In order to more closely examine sources of instability in the 
above type of problem, and to determine the proper variable for a root 
locus study, it is useful to solve for the system* s uncoupled form |~3 
Rewriting ( 3 ) in the time domain gives 



h(t) + 2 ah(t) + w 2 h(t) = KVd(t) 



(8) 



li 



In order to rewrite (8) in matrix form let 

x 1 (t) = h(t) 

= h(t) 

• M 

= h(t) 
x^(t) = h”(t ) 
u(t) = KVd(t) 

This leads to the following matrix form for the system; 



( 9 ) 



Xg(t) 

* 3 (t) 





0 


l 


0 


0 




0 


• 


0 


0 


1 


0 




0 


X = 


0 


0 


°p 


1 


X + 


0 




0 


0 


2 

-w 


-2a 




1 



( 10 ) 



which is defined as 

x = Fx + Du ( 11 ) 

The system can now be transformed into the uncoupled form by 



defining a new variable y as 
y L = u(s)/s 2 
y 2 = u(s)/s 

2 2 
= u(s)/(s + 2as + w ) 

2 2 

y^ = u(s)s/(s + 2as + w ) 

After expanding x by partial fraction expansion it is noted that 



( 12 ) 



x = 



l/w^ 

0 

0 

0 



•2a/w 2 

1/w 

0 

0 



(4a 2 -w 2 )/w2 

-2a/w 

1 

0 



2 a /wp 

- l/w 
0 
1 



( 13 ) 



which is defined as 
x = 



( 14 ) 
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Also, by solving for the inverse of G, the expression 



-1 

^ = G x 

can be written as 



1 = 



w 

0 

0 

0 



2 2 

w 

0 

0 



1 

2a 

1 

0 



0 

1 

0 

1 



( 15 ) 



( 16 ) 



In order to draw a signal flow diagram of the uncoupled system,, 
(n.) is substituted into (ll) giving 
Gj£ = FG^ + Du 



(17) 



and multiplying by G 

„-l__ -1 

^ = G FG^ + G Du 

Solving (18) gives the uncoupled form 

0 

1 
0 

1 





0 


1 


0 


0 




* 


0 


0 


0 


0 




1 = 


0 


0 


Op 


1 


1 + 




0 


0 


-W 2 


-2a 





(18) 



(19) 



which is drawn as shown in Figure 2. Lines are added to denote the 
subsequent addition of state variable feedback loops. 
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